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1 Abstract 

Thermal neutrons are of considerable interest to the Department of Defense and for 
commercial applications. Unlike high-energy photons, neutrons easily penetrate high 
density targets, but get effectively absorbed by low density materials like paraffin, nylon 
or explosives. This makes them attractive complements to X-rays for radiographic 
applications, e.g. for the detection or inspection of explosives inside steel casings. The 
key challenge is to develop a compact generator for thermal neutrons with large enough 
flux. The limited availability of radio-isotopes, combined with the relatively short half- 
life, safety constraints and regulatory requirements make them unattractive for wide¬ 
spread use. An alternative design exploits the Deuterium-Tritium (D-T) fusion, which 
generates Alpha particles and fast neutrons. In these sources, Deuterium ions are 
accelerated to about 130 keV and hit a Tritium target. The acceleration of Deuterium ions 
is usually accomplished in a diode configuration. Recently, considerable success has been 
achieved in the acceleration of ions via laser-matter interaction. In this project we 
investigated whether laser-accelerated ions could undergo nuclear fusion in an adequately 
designed target and could be used for neutron sources. We therefore enhanced our 
proprietary plasma simulation code VORPAL with a model for fusion reactions and 
investigated the generation of neutrons in shaped D-T targets. We find that neutron fluxes 
large enough for radiographic applications can be generated by utilizing moderate (~ 
10 17 —10 18 W/cm 2 ) laser intensities. 
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3 Introduction 

Thermal neutrons with energies on the order of 0.025eV are of considerable interest to 
the Department of Defense and for commercial applications. Unlike high-energy photons, 
neutrons easily penetrate high density targets, but get effectively absorbed by low density 
materials like paraffin, nylon or explosives. This makes them attractive complements to 
X-rays for radiographic applications, e.g. for the detection or inspection of explosives 
inside steel casings. The key challenge is to develop a compact generator for thermal 
neutrons with a flux of at least 10 8 neutrons/cm 2 s which is sufficient for radiographic 
applications. 

The limited available amount of radio-isotopes like Californium 252, combined with the 
relatively short half-life, safety constraints and regulatory requirements make them 
unattractive for wide-spread use. 

An alternative design exploits the Deuterium-Tritium (D-T) fusion, which generates 
Alpha particles and fast neutrons. In these sources, Deuterium ions are accelerated to 
about 130 keV and hit a Tritium target. In order to be attractive for radiographic 
applications, the 14.1 MeV neutrons have to be thermalized in an external moderator, for 
example a paraffin layer. This reduces the useful neutron flux by about two orders of 
magnitude. D-D reactions produce neutrons with lower energy (2.4 MeV), but one must 
also account for the branching ratio into charged products, i.e. d(d,p)t, and the reduced 
cross-section compared to d(d,n)a. 

The acceleration of Deuterium ions is usually accomplished in a diode configuration. 
Recently, the interaction of high-power laser pulses with plasmas or solids has attracted 
considerable interest for the acceleration of electrons or ions [1, 2], Especially the 
acceleration of ions to MeV energies opens entirely new possibilities for small scale ion 
accelerators and possibly for neutron sources. 

In this project, we therefore investigated different target configurations for the potential 
use as neutron generators via numerical simulations. 


4 The plasma simulation framework VORPAL 

The versatile plasma simulation code VORPAL [3] was originally designed for the 
investigation of laser wakefield acceleration of electrons. However, the code design is 
sufficiently general that it is capable to address entirely new application areas, like 
modeling of superconducting RF cavities, breakdown phenomena in microwaves guides 
or investigations of astrophysical plasmas. Figure 1 shows examples of such simulations. 
The code has recently been enhanced by a variety of models, including ionization 
processes. 
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The model for ionization processes is a hybrid between classical particle-in-cell [4] and 
the Direct-Simulation Monte-Carlo model [5]. Particles are pushed in their self-consistent 
electromagnetic field. At every time-step, particles within a cell are considered for 
possible collision (and therefore reaction) by Monte-Carlo methods. The main parameters 
required in this model are energy dependent cross sections for ionization. While the 
model was mainly designed for ionization processes, it can relatively easily be adapted to 
other reaction types, including fusion. 



Figure 1: Left panel: Example of a VORPAL simulation, (superconducting RF cavity). Right Panel: Mono- 
energetic electron beam generated by laser-wakefield acceleration (Simulation by John R. Cary, 
Visualization: Peter Messmer). 


4.1 Ionization Model 

The ionization model in VORPAL was successfully used to reproduce ion spectra in 
Electron-Cyclotron Resonance Sources. Figure 2 shows a comparison of the kinetic 
impact ionization model in VORPAL with a spectrum for a hot population of oxygen 
measured at VENUS [7], For experimental reasons, the spectrum was not measured for 
pure oxygen, but rather a mixture of Oxygen and Helium. The resulting spectrum at 
mass/charge ratios of 4 was therefore a blend between the oxygen and the helium 
contribution. The VORPAL simulation is a fully periodic box, containing atomic 
Oxygen, 0 + and e". The evolution of the different charge states are shown in Figure 4 
(right). For the given electron energy, the first charge state not present in the initial 
distribution is going to be 0 2+ , but as soon as it is available in high enough concentration, 
the next ionization state is being produced and the ionization rate of 0 2+ drops. The 
overall spectrum obtained by integrating the ionization rates in time is shown in the left 
panel. For comparison, the measured current at VENUS [8] is shown. Good agreement is 
obtained for 5 charge states. As mentioned previously, the charge state M/Q=4 is mixed 
in the measured spectrum. We therefore plotted a token value of 0 at M/Q=4. The good 
agreement of the simulated with the measured spectrum indicates that the kinetic 
ionization model in VORPAL works as expected. 

The energy dependent cross-sections for the ionization model in VORPAL are based on 
the semi-empirical model by Shull and Van Steenberg [9], which originates from the 
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astrophysics community. This models predicts the energy dependent cross-sections of 
light elements, including, C, N, O, Ne, Mg, Si, S, Ar, Ca, Fe, and Ni. 



Figure 2: Left: Comparison of an oxygen spectrum computed with the kinetic ionization model in 
VORPAL [bars] with an observed spectrum [red stars]. The measurement at Oxygen 4+ was blended with 
Helium and therefore not shown in the picture. Right: Temporal evolution of the different charge states of 
Oxygen for a population starting with pure Oxygen, and 0 + : 0 2+ (black), 0 3+ (brown), 0 4+ (organge), 0 5+ 
(yellow), 0 6+ (purple), 0 7+ (blue). 


4.2 Fusion Reaction 


The first part of this project was to implement a cross-section model for fusion reactions 
into VORPAL. Due to the relatively large cross-section at low energies, we chose the 
Deuterium-Tritium (D-T) reaction for all further studies. We chose a parameterized 
cross-section model [6] of the form 


v T (E) = 


A 5 + 

(A 4 -A 3 E) 2 + if 


E 

exp {A,E- V2 )-\ 

r 


with Ai = 45.95 barn keV, A 2 = 50200 barn keV, A 3 = 3.98e-3 keV' 1 , A 4 =1.297 and 

-24 2 

A 5 =409 bam keV, and the total cross section Oj measured in barn (10‘ cm ) and the 
energy E in keV. 

Figure 3 shows the cross-sections as a function of energy. The peak cross-section is on 
the order of 10~ 27 m 2 , centered around an energy of 130 keV. 

One of the main differences between the kinetic ionization model and the fusion model is 
the magnitude of the cross-section, which is about 7 orders of magnitude smaller for 
fusion than for ionization. We therefore did not track the generated 4 He (a particles) in 
the simulation, but rather estimated the number of neutrons generated in the interaction 
based on the particle statistics. This de-coupled approximation is valid because the a- 
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particles are too few to have an impact on the plasma dynamics (i.e. energy deposition), 
despite the short mean free path at such high densities. 


4.3 Fusion of shock-accelerated ions 

It has previously been reported that the interaction of a femto-second laser pulse with 
over-dense plasmas can lead to the formation of a shock wave in the target and to the 
acceleration of ions at this shock front [10]. In a target consisting of a D-T mix, this will 
therefore lead to a population of D-T ions streaming through a D-T background. If the 
energy of the shock- accelerated ions is large enough compared to the background ions, 
one expects that fusion reactions and therefore neutron generation will occur in the 
vicinity of the shock front 

In order to test this hypothesis, we modeled the interaction of a femto-second laser pulse 
with a slightly over-dense D-T target. We assumed a plasma of density 2.8 10 27 m' 3 being 
irradiated by a 60 fs, 0.8 pm laser pulse with an intensity of 1=1.8 10 17 W/cm 2 and a spot 
radius of 3 pm. The plasma is over-dense with a ratio of n/n cr = 1.6. The plasma was 
modeled in a 2D configuration on a 3000 x 5 cell grid. The grid resolution in longitudinal 
direction is 1.6 nm and 10 nm in transverse direction. Perfectly matched layer (PML) 
boundary conditions were employed along both directions of propagation and periodicity 
was assumed in the transverse direction. 


D-T Crossections 



Figure 3: Energy dependent cross-section for the D-T fusion reaction based on the model in [6]. The peak 
cross-section is around 100 keV and reaches 10' 27 m 2 (10 bam). 


Figure 4 shows phase space projections of the interaction at different times. Due to the 
lighter mass, the Deuterium ions get accelerated to higher velocities than the Tritium 
ions. The Target Normal Sheath Acceleration (TNSA) mechanism accelerates both D and 
T ions at the backside of the plasma (at x=2.5 pm). However, in addition to the TNSA, 
particles are accelerated in the interior of the plasma at the shock front. This results in a 
fast population of D (and T) propagating forward through the background plasma. 
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Figure 4: Phase space projections of the D (blue) and T (red) ions at different times: (a)- 96fs, (b)-0.48ps 
and (c)- 0.96ps. 
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Figure 5 shows the energy distribution of the ions for the above parameters. While these 
simulations demonstrate that a fast population of D can be generated by the laser-plasma 
interaction, the resulting ion energies are only on the order of 50 keV, for which the 
cross-section (see Figure 3) is too small to generate any significant number of fusion 
events. 




Figure 5: Spatial energy distribution of the D (blue) and T (red) population after 0.48ps (left) and 0.96ps 
(right) for a laser intensity of 1=1.8 10 17 W/cm 2 . 


Figure 6 shows the ion phase space and energy of an interaction with increased laser 
intensity of I = 2.6 10 17 W/cm 2 . All other parameters are identical to the previous 
simulations. The energies of the shock-accelerated ions is well above the 120 keV 
required for a significant number of fusion events to occur. 



0.5 1.0 1.5 2.0 2.5 3.0 0.5 1.0 1.5 2.0 2.5 3.0 

x [/^m] x [j&nn] 


Figure 6: Phase space projection (left) and ion energy (right) for the interaction of a 60 fs, 2.6 10 17 W/cm2 
laser pulse with a D-T plasma. 

Using the newly developed fusion model in VORPAL, we were able to investigate the 
neutrons generated in this interaction. Figure 7 shows the temporal evolution of these 
neutrons. The peak neutron flux occurs when the shock-accelerated ions propagate 
through the unperturbed plasma. Once they reach the former back side of the plasma, 
which is now accelerated via TNSA, the relative velocity is reduced and the number of 
fusion events is therefore reduced. 
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Figure 7: Left: Temporal evolution of the neutron flux generated in the interaction of a 60 fs laser pulse 
with a D-T target. Right: Scaling of the peak neutron flux with increasing laser intensity. 

The energy of the shock-accelerated ions is highly dependent on the shock velocity, 
which is in turn dependent on the laser intensity. As these ions are responsible for the 
neutron generation, the laser intensity will control the neutron flux. Figure 7 shows the 
scaling of the peak neutron flux as a function of laser intensity. At intensities > 3 10 17 
W/cm 2 , the target was destroyed. 

In summary, using the newly developed fusion model in VORPAL, we demonstrated that 
the interaction of high-intensity laser pulse with an over-dense D-T target can lead to the 
generation of neutrons. 



- 1.5 

- 2.0 


0 10 20 30 40 50 

X [/Linn] 

Figure 8: Reverse acceleration of ions in a two material configuration. Laser pulse entered from the left 
and accelerated the first population. At the interface of the two materials, plasma gets accelerated in 
backward direction. 
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4.4 Reverse acceleration of ions 

While we demonstrated in the previous paragraph that the neutron flux generated by the 
interaction of a femto-second laser pulse with a D-T mix is sufficient for radiographic 
applications, one of the drawbacks of that configuration is the significant fraction of the 
D-T mix accelerated in the reverse direction (i.e. opposite to the initial laser propagation 
direction) which is lost to neutron generation. These back-propagating particles can be 
seen in the phase-space plots of Figure 4b-c, for x < 1.5 pm. A more efficient target 
design would take advantage of that ablated material as well. 

One way of taking advantage of the ablated material is to optimize the target such that 
ions are only accelerated in the reverse direction. This can be accomplished by a target 
consisting of an under-dense coating applied to an over-dense material. 

Figure 8 shows the ion phase space of such a two-component target, demonstrating the 
ion acceleration mainly in the reverse direction. This is an interesting and potentially 
important mechanism for this type of target, which was not previously investigated. 

This mechanism is not limited to pseudo ID simulations, but also works in 2D. Figure 9 
shows the propagation of a 10 21 W/m 2 , 60 fs, 0.8 pm laser pulse in a Helium coated 
aluminum target. The simulation was run on a 640 x 320 cell grid, spanning a 64 x 32 pm 
domain. The neutral Helium is modeled as a gas at density of 10 26 #/m" 3 . The propagation 
of the laser pulse ionizes the coating and eventually accelerates the ions. This 
demonstrates that the reverse acceleration of ions also works in higher dimensions. 



Figure 9: Ion phase space, x-y-v x , showing an iso-density surface for He 2+ ions (yellow) and Al + ions 
(green). Negative ion velocities are plotted upwards. The laser pulse initially propagates in positive x 
direction, gets reflected off the A1 mirror and leaves the domain in negative x direction. The He 2+ ions 
initially coat the inside of the spherical A1 mirror and are accelerated in reverse direction by the laser pulse. 
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4.5 Focusing of reverse accelerated ions 

In this section we therefore investigate whether the material ablated in the interaction of a 
strong laser pulse with an over-dense target can lead to a neutron flux sufficiently large 
for radiographic applications. An alternative application for these reverse-accelerated 
ions would be to use them as a driver for e.g. Inertial Confinement Fusion. In particular, 
one may envision a combination of forward-accelerated (sheath mechanism) and reverse- 
accelerated ion beams created by the same laser pulse, operating on different sides of a 
target, to obtain synchronized beam-beam interactions. Additionally, one may also 
envision the target as spherically shaped, resulting in converging beams to increase 
density and temperature at the beam focal point; the target would then include holes 
through which the laser beam can pass, similar to Hohlraum targets. 

We therefore start with investigating the possibility of converging ion beams from shaped 
targets. We ran simulations of a 10 pm thick foil, with a "dimple" on the near side having 
a radius of curvature of 5 jam. The plate is comprised of a neutral plasma with equal 
numbers of electrons and singly-ionized deuterium ions. Both plasmas have a density of 

2.5 10 27 m' 3 . The plasma frequency is C0p=2.825 10 15 rad/s and the initial temperature is 
130’000 K, resulting in a Debye length of 0.866 nm. While this leads to numerical 
heating of the plasma, we anticipated that the violent laser-plasma interaction will 
dominate and provide at least some qualitative results. For more quantitative results, we 
then performed subsequent simulations at higher resolution. 



6xi <r B 

4x10" 0 

2x10" B 

0 

— 2x1 Q -0 

- 4x10 -B 

-6X1G" 0 

4.QX1O -0 S-QxlO" 0 0.0*1 Q -6 1.0x10" 5 1.2*1 Q -5 l,4xlCT 5 1.0x10” 


Figure 10: Flow field of reverse accelerated ions. 


The laser has a wavelength of 0.8 pm, for which the plasma is slightly over-dense, 
n/n cr it=1.44. The laser field has a peak amplitude of 1.5T0 12 V/m (corresponding to an 
intensity of 3T0 17 W/cm 2 ) and a pulse duration of 20fs. 

These initial simulations were run on a 3000x2000 cell grid spanning a physical domain 
of 30 pm x 20 pm (dx = dy = 10 nm), and a time-step of 1.73 10" 17 s. Perfectly matched 
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layer boundary conditions were used along the laser propagation direction and periodic 
boundaries in transverse direction. 

Figure 10 shows the Deuterium flow field after 1.73 ps (lOO’OOO timesteps). The figure 
clearly shows the backward motion, as well as a focusing flow of the ablated material, 
which could eventually lead to a high density population in the center of the ‘dimple’. 

One of the key problems encountered in these simulations is the long timescale on which 
the interaction takes place. In addition, it would be highly useful to be able to fully 
resolve high density plasmas to avoid numerical heating, a common problem in Particle- 
In-Cell (PIC) simulations when the Debye length is not well resolved. In order to further 
investigate the focusing and to actually observe the formation of the focused ion 
population, we ran additional simulations, now with a finer grid resolution of 6000x4000 
cells (dx =dy = 5nm) and a time-step of 1.17 10~ 17 s. In addition to the finer grid 
resolution, we also applied a [1 2 1] box-car filter to smooth each component of the 
current deposited by the charges at every time step. This resulted in a significant 
reduction of numerical heating and therefore enabled simulations on long timescales. 
Finally, we increased the laser pulse duration from 20 fs to 50 fs. 



Figure 11: Two different views of the ion density after 1.18 ps. The laser pulse entered from the left, was 
reflected at the spherical mirror and left the domain to the left. The reversely accelerated ions form a high 
density population in the center of the mirror. The density of this population is significantly higher than the 
TNSA accelerated population on the right hand side of the foil. The boundaries of the initial distribution are 
shown as dotted lines. 


The basic steps of the laser-matter interaction in these fully 2D simulations are essentially 
the same as in the pseudo ID simulation in section 4.3, i.e. TNSA. In this mechanism, the 
electric field of the laser pulse accelerates the electrons in the interaction region, which 
are rapidly accelerated into the target by the Lorentz force. The electrons stream through 
the foil and form a high-temperature, charged cloud out on the back side. The ambipolar 
field created by this cloud (sheath) accelerates the ions forward, while colder electrons 
are pulled back into the material (return current). If the laser pulse is long enough, the 
electrons return into the interaction region while the laser pulse is still on, being injected 
again into the foil (electron recirculation). At the same time, the ions in the interaction 
region start to move in reverse direction, due to the local positive charge left inside the 
target (reverse accelerated ions). 
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While ID simulations capture all these effects, they are unable to answer questions about 
the spatial distribution of the reversely accelerated ions, especially in a shaped target. 
They also do not properly capture effects from transverse gradients; for example, the 
return current at the periphery of the channel (outside the laser pulse envelope) can 
generate strong magnetic fields and affect the ion propagation. 



Figure 12: Evolution of the ion density over time. The laser pulse propagated along the x direction and got 
reflected at the spherical mirror. Time evolves in negative Y direction. The ‘ripples’ on the ‘canyon-wall’ 
show the strong filamentation of the ablated material. After about 1 ps, a high density population in the 
center of the mirror has formed, propagating in backward direction. This population shows strong 
divergence, most likely due to the high local space charge. 


4.6 D-T acceleration in Shaped Targets 

While strong filamentation of the reverse-accelerated ion beam may prevent a good 
focusing for beam-beam interactions, we still need to examine the potential use of the 
reverse-accelerated material shown in the previous section as a driver for nuclear fusion. 
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We test a generic configuration where the reverse-accelerated ions can be focused onto a 
target. Figure 13 shows the overall setup: the main target consists of a 10 pm thick tin foil 
at a density of 2.5 10 27 m' 3 with a 5 [tin hemi-spherical cavity in the center, which acts as 
a mirror. The cavity is coated with an under-dense, 2.5 [tin thick Deuterium layer at a 
density of 4.4 10 26 m~ 3 . Suspended at the focus of the mirror is a Tritium ball with 2 jam 
diameter, also at a density of 4.4'10 26 m~ 3 . 

The 0.8 jam laser has a spot radius of 2.5 jam and an intensity ranging between 1.5 10 21 
Wm' 2 and 2.4 10 22 Wm' 2 corresponding to peak electric fields of 1.5 TV/m to 6 TV/m. 
The simulation domain is 30 x 20 jam covered by a computational grid of 6000 x 4000 
cells and a simulation time-step of 1.0610~ 17 s. Perfectly matched layer boundary 
conditions are applied in the laser propagation direction and periodic boundary conditions 
in the transverse direction. Similar to the simulations in the previous section, we chose an 
electron temperature which results in a Debye length resolved by the computational grid, 
as well as filtering of the currents in order to suppress numerical noise. 



Figure 13: Setup of shaped target configuration: A high density shaped tin target (red) is coated with an 
under-dense deuterium film (yellow) with a Tritium ball (blue) suspended in the center of spherical mirror. 
The Tritium ball has a diameter of 2 pm, the Deuterium coating has a thickness of 2,5 pm. The 
cavity/mirror radius is 5 pm. The laser pulse enters from the left, is reflected off the spherical mirror and 
gets focused in the center of the mirror. 


Figure 14 shows the development of the high density jet of ablated Deuterium. The over- 
dense lead foil, coated with the deuterium, is not shown. The peak density in the ejected 
Deuterium can reach up to 1.5 times the original Deuterium density. The Tritium ball 
explodes as a result of expulsion of electrons by the laser pulse (Coulomb explosion) and 
interacts with the counter-streaming Deuterium. 
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Figure 14: Time evolution of the D and T density for a laser amplitude of E=3el2V/m (left) and 6el2 V/m 
(right) after O.lps (a), 0.3ps (b), 0.6ps (c) and 0.9ps (d). 
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Figure 15: Deuterium (left) and Tritium (right) density distribution after 1.33ps for a peak electric field of 
E= 3TV/m. 


Figure 15 shows the Deuterium and Tritium density distribution after 1.33ps. While at 
that time the reverse accelerated ions form a high density peak at the original location of 
the Tritium ball, the Tritium has been spread over a large volume due to Coulomb 
explosion. 


4.7 Neutron Generation in Shaped Target 

While the previous paragraph has shown that the ablated material in a shaped target can 
be focused and generate a high density peak, it remains to be seen whether the 
accelerated ions have the right energy to undergo fusion. 

Figure 16 shows the temporal development of the neutron generation inferred from the 
D/T densities and local mean velocities for two different laser amplitudes. For both 
amplitudes, neutron generation initially happens where the ions from the Coulomb 
explosion of the Tritium ball penetrate the Deuterium. Later, the interaction of the 
Deuterium jets creates another location of high neutron yield. Finally, the Tritium density 
at the original Tritium ball location becomes so low, that neutrons are no longer produced 
there. 

Integrating the number of generated neutrons over the entire simulation domain allows 
estimating the total neutron flux as a function of time. Figure 17 shows the temporal 
evolution of the neutron flux for different laser peak electric fields. While increasing the 
laser peak field leads quicker to the generation of neutrons, the ions exceed the drift 
velocity of peak fusion cross-section. They therefore are lost to the generation of 
neutrons. On the other hand, for too weak a laser field, the accelerated ions are not fast 
enough and therefore do not lead to nuclear fusion. For the given configuration, a peak 
electric field of 3 TV/m maximizes the neutron flux. 
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Figure 16: Time evolution of the neutron generation for a laser amplitude of E=3el2V/m (left) and 6el2 
V/m (right) after 0.3ps (a), 0.6ps (b) and 0.9ps (c). 
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Figure 17: Neutron flux evolution for laser amplitudes E=1.5 TV/m (dash-dotted), E=2 TV/m (dash-triple- 
dotted), E = 3 TV/m (solid) and E = 6 TV/m (dashed). Increasing the amplitude leads to an increased 
momentum of the accelerated ions and therefore to a relative ion energy which is beyond the peak cross- 
section for the D-T fusion reaction. 


The neutron flux in Figure 17 is computed at a distance of 1 m away from the target. 
Even assuming a reduction of the flux by about two orders of magnitude due to an 
external moderator should lead to a sufficient flux of neutrons for radiographic 
applications. 

4.8 Summary 

In summary we found that ions ablated in the interaction of a strong laser pulse with an 
over-dense plasma can be focused by shaping the target. In case of a spherically shaped 
target, this leads to the formation of a high density ion population in the center of the 
spherical mirror. Target shaping has been successfully used to focus ions accelerated via 
TNSA on the backside of the target. These accelerated ions are then used as drivers for 
Inertial Confinement Fusion. Our simulations indicate that a similar population can be 
generated at the interaction site, leading to a drive beam traveling in the backward 
direction. 

We have shown that the interaction of a femto-second laser pulse with a slightly over- 
dense target consisting of a D-T mix can be used to generate a neutron flux sufficiently 
large for radiographic applications. However, a large fraction of the D-T mix does not 
participate in the D-T fusion reaction and is therefore lost. 

In a last configuration, the laser pulse was used to strip a Tritium ball from all its 
electrons. The resulting Coulomb explosion was used to accelerate the Tritium ions. The 
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reverse ion acceleration mechanism in a shaped target generates a counter-streaming 
population of Deuterium ions. Depending on the peak electric field of the laser pulse, a 
neutron flux large enough for radiographic applications can be generated. 

The above simulations pushed 40 million particles on a 6000x4000 cell grid for about 
110’000 time-steps and produced about 250 GB of data and took about 500 CPU days 
(12’000 CPU hours). Simulations of this size are too large for performing e.g. Parameter 
studies. One possibility to reduce the simulation time is to employ higher-order particle 
shapes, which interpolate the electromagnetic field not only from the surrounding cell, 
but from a more extended region in space. Using the field values from three neighboring 
cells, it is possible to relax the constraint to resolve the Debye length and only requiring 
to resolve the electron inertial length. This results in a grid resolution which is lower by a 
factor (c/vtherm) 3 which can easily exceed a factor of l’OOO. In addition, this yields to an 
increased maximum possible time step due to its coupling via the CFL condition. Time 
savings of factors 10’000 and more can therefore be achieved. 


5 Conclusions 

In this project, we have implemented a kinetic model for nuclear fusion in the plasmas 
simulation code VORPAL. We therefore extended a kinetic reaction model implemented 
in a different project by cross-sections for the Deuterium-Tritium nuclear reaction. 

We investigated the interaction of a medium powered laser pulse with a slightly over- 
dense thin Deuterium/Tritium target. The interaction leads to the formation of a shock 
wave inside the plasma, which can accelerate particles at the shock front. This leads to 
the formation of a population of D (and T) streaming relative to a D/T background. 

Using our newly developed fusion reaction model, we were able to compute the neutron 
flux generated in this interaction. The peak neutron flux was on the order of 10 12 #/cm 2 s, 
with a strong dependence on the incident laser pulse intensity. 

While the overall number of particles per shot is relatively small due to the short 
interaction time, the high repetition rate of commercially available lasers with the 
required parameters can lead to high average neutron numbers. 

Discussions with experimentalists have indicated that slightly over-dense thin targets can 
be realized using foams or gels or even using high-density gas jets. The simulation results 
therefore could be tested experimentally and possibly lead to the development of a 
compact neutron source. 

We also performed simulations of reverse acceleration by laser-plasma interaction in 
shaped targets. We have demonstrated that a focusing of the ablated material occurs and 
leads to a high-density ion population propagating in reverse direction. This population 
could be used as a driver for Inertial Confinement Fusion. 
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Particle-in-Cell simulations of high-density plasmas are computationally highly 
demanding and algorithmic enhancements should be investigated. For example, higher 
order particle shapes have demonstrated that considerable time savings can be achieved 
compared to first order interpolation for particles. 
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